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I investigate stable equilibria of bodies in potential fields satisfying a generalized Poisson equation 
V • [p(\ V<p\/ao)V<p] oc p. This describes diverse systems such as nonlinear dielectrics, certain flow 
problems, magnets, and superconductors in nonlinear magnetic media; equilibria of forced soap films; 
and equilibria in certain nonlinear field theories such as Born-Infeld electromagnetism. Earnshaw's 
theorem, totally barring stable equilibria in the linear case, breaks down. While it is still impossible 
to suspend a test, point charge or dipole, one can suspend point bodies of finite charge, or extended 
test-charge bodies. I examine circumstances under which this can be done, using limits and special 
cases. 1 also consider the analogue of magnetic trapping of neutral (dipolar) particles. 
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$_i ■ I. INTRODUCTION 



The possibility to suspend bodies, or to levitate them against the pull of gravity, by applying various long-range 
fields to them, is of obvious interest and importance. This has to be done by fields whose sources are outside the 
body to be suspended. The fields usually considered (electrostatic, magnetostatic, nonrelativistic gravitational) are 
derived from scalar potentials that satisfy the Laplace equation outside sources. There is a sweeping statement, under 
the name of Earnshaw's theorem [jj], of the impossibility to suspend, statically, in stable equilibrium, a body carrying 
any combination of the different charges (in rigid distributions) in any setting of the fields, when there is no overlap 
between suspended charges and sources. This has to do with the fact that a solution of the Laplace equation can 
^C) i attain true extrema only on the boundary of the domain of solution-the maximum (extremum) principle, see e.g. [Q. 

Many of the Earnshaw-type results, beyond the basic statement of the impossibility to suspend a point charge, rest 
heavily on the linearity of the Laplace operator. The Earnshaw statement may thus be expected to break down in 
non-linear systems. Here I revisit the question of stable equilibria of bodies subject to fields described by a nonlinear 
generalization of the Poisson equation: A source distribution p(r), in D-dimensional Euclidean space, produces a 
potential if through 
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■ A/V = V • [n(\V<p\/ao)$<p] = a D G P , (1) 

s : 

O ■ derivable from the action 

o ■ 2 

S = -Jppd°r- ^ J T[{V^f/a 2 ] d°r. (2) 

'X: - 

Here, ao is a constant with the dimensions of Vtp, G is a coupling constant, a D — 2(n) D/2 /T(D/2) is the D-dimensional 
solid angle (introduced for convenience), p(x) = dJ-(y)/dy {y — x 2 ) is positive except perhaps that p(0) may vanish. I 
normalize to ^"(0) = 0. The nonlinear operator Af generalizes the Laplacian, and can also be written as Af — pAijdidj , 
with Aij = S{j + pap^ipj /\Vip\ 2 , where p = xp! j p. Points where Vip = need special treatment in theories where 
p(0) = 0, but this does not modify any of the results, and is ignored for brevity. (Summation over repeated indices 
is implied throughout.) Only the case where eq.(Q) is elliptic is considered, so Aij is positive definite, tantamount to 
fx > — 1. For G > 0, like, point sources attract each other, and opposite sources repel (as in gravity), and vice versa 
for G < (as in electostatics) H. 

In H I derive general results pertaining to forces on bodies in such theories. Here I concentrate on the existence 
of, and criteria for, stable equilibria of bodies. Solving the general problem requires numerical means, but much can 
be learnt from closed-form solutions that can be found in certain limits and special cases. 

Equation(Q) describes a variety of physical problems. Some examples are: (a) Nonlinear dielectric, and diamagnetic, 
media; p is then the dielectric, or diamagnetic coefficient, which depends on the field strength (G < 0). (b) Subsonic- 
potential- flow problems of non- viscous fluids with an equation of state of the formp = p(g) (subsonicity is equivalent to 
ellipticity). The form of p(x) depends on the equation of state. For instance, when p cx g 7 , p(x) on [1 — (7 — l)x 2 ] 1 ^ 7-1 ^ 
for 7 > 1 (the limit 7 — > 1 exists: exp(— x 2 )) (G > 0). Our results then pertain to equilibria of sources, sinks, and 
obstacles in the flow, (c) Equation^) was used in Q| as an effective-action approximation to Abelianized QCD, with 
p(x) oc Inx; it is elliptic for x < e" 1 , or x > 1. (d) Nonlinear (vacuum) electrostatics as formulated e.g. in the 
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Born-Infeld nonlinear electromagnetism, which also appears in effective Lagrangians resulting from string theory (see 
review and references in |J jfj)). In the electrostatic case p{x) oc (1 — x 2 ) -1 / 2 , and G < 0. (e) A formulation of an 
alternative nonrelativistic gravity to replace the dark-matter hypothesis in galactic systems 0. Here p,(x) ~ x for 
x<l, and \i « 1 for i > 1 (G > 0). (f) Problems of nonlinear electric-current flows in systems with field-dependent 
conductivity (nonlinear current-voltage relation), and nonlinear diffusion problems; p{x) is the transport coefficient. 
Here, a force on a body signifies the gradient of the entropy-production rate with position of the body, so stable 
equilibria are configurations of extremal entropy generation rates, (g) Area (volume) minimization problems: If p(r) 
is understood as the height of a D-dimensional surface above position r on a Z?-dimensional hyperplane H, then eq. (jl|) 
describes the problem of the minimization of the volume of the surface. The sources could represent a prescribed 
vertical- force density. Forces on sources as will appear below correspond to lateral forces (parallel to H). In this 
problem p,(x) oc (1 + a; 2 ) -1 / 2 , and G > 0. 



II. THE GENERAL PROBLEM 



I start with the analogue of Earnshaw's original question: can an equilibrium configuration p(r) exsit subject to the 
tp field alone? This would require that pVp vanish everywhere, and can probably be precluded, as I show for a wide 
class of theories. In || I derive an expression for the virial integral V = f pr ■ Vp d D r — a^/2a D G J J-(D — 2T) dPr 
[with J-(y) = yJ- 1 (y) / J-(y)], that holds for theories in which J-(0) < D/2. (In this case the potential vanishes 
asymptotically for a bounded charge; for J-(0) = D/2 it diverges logarithmically.) We learn from this that in 
problems with T{y) < D/2 we have V ^ 0, and /cV(/j cannot vanish everywhere (actually true even for T = D/2 
where the expression for V has an extra term). This applies e.g., to the flow, and volume- minimization, problem 
(where T < 1), and to the modified dynamics (where 1 < T < D/2). Such complete-equilibrium configurations 
seem, anyway, to be of academic interest only. The bodies making up the systems we want to keep in equilibrium are 
themselves not held together, internally, in static equilibrium by the same forces: atoms and stars consist of moving 
constituents, held, intrinsically by forces other than electricity, and gravity, respectively. 

Consider then equilibria of a charged body B, given by a rigid charge distribution p B (r), in the presence of some 
fixed distribution p(r) of "holding" sources, with no overlap: pp B — 0. (Body and sources might be each held rigid 
by forces other than the p field.) The force on B is 

p B Vip d D r, (3) 

with p determined from eq.(|l]) with p(r) — p B (r) + p(r). F is also the gradient of the energy with respect to rigid 
translations of B by 5a: 5E = — 5a- F [3]. I also want to consider the possibility of several types of charge p a coupled 
to fields p a through equations of type (1), possibly with different forms of p{x) (e.g. levitation in some tp field against 
gravity). Under p — > — p, p — > — tp , and F->F; under G — > — G, tp and all forces change sign. 

Under what conditions (choice of p(x), boundary conditions, p, and p B ) can B be suspended stably (with £ at a 
minimum)? I am mainly concerned with translational stability, so define -E(R) as the translational enery as a function 
of the position of some reference point R in B, keeping the orientation fixed. F = — Vr-E(R). A stable equilibrium 
where, F = 0, requires that the Hessian of E, = d ^.g R . , be positive definite. (More generally, that the first 
non-vanishing derivative-which must be of even order-is positive definite: E :il .. %i2 aj 1 ...aj to > for any non-zero a.) 

In the linear case we can separate tp — tp + p B , with tp, and p B coming, respectively from the "holding" sources, 
and from B alone. The contribution of p B to the force in eq.(^) vanishes, because a body does not exert a force on 
itself, and we then have (for the many-field case) 

E(B.)=J2 /^«(R + x)(to (4) 

a J 

(x is the position inside the body). Clearly, when all tp a satisfy the Laplace equation, so does E(R), and E does not 
have a true extremum, leading to the general Earnshaw statement. 

In the nonlinear case eq.(Q) remains valid only when B is a test body, whose contribution to the sources can be 
neglected when calculating the potential. For a point, test body, in a single field we have -E'(R) = qp>, and it is still 
not possible to have a stable equilibrium because solutions tp of eq.(Q) still satisfy an extremum principle. However, 
any departure from this restricted case might permit stable suspension, (i) Unlike a test point charge, a point body 
of finite charge is not denied a stable equilibrium. The extremum principle for tp does not exclude this, as tp is no 
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more the effective potential (energy) of the body, (ii) Even test charges arranged rigidly in an extended body can 
sometimes be suspended in a tp field, (iii) A point, test particle that carries test charges of different types can be 
suspended in a combination of fields tp a satisfying eq. ([!]). Now, E = J2 a Qa^a is not subject to an extremum principle 
even though the separate tp a s are. 



III. EQUILIBRIA OF POINT CHARGES 

Take the case of a point charge q "held" by a distribution p, when the only boundary condition dictated is tp —> 
at infinity. The problem can be solved in closed form in the limit where q is very large (relative to p). The total force 
on the whole system vanishes, so the force on q is equal and opposite that on p. In the above limit, p may be treated 
as a distribution of test charges, so the force on it is — J d D rp(r)Vip B , where (p B is produced by B alone. This is 
gotten straightforwardly by applying Gauss theorem to eq. (jl|) for a point charge. Putting all this together we obtain 
for the force F g (R) on q, at position R: 

P,(R) = s( qG )a J d °rp(r)S ( J^ ) ^ (5) 

Here, s(x) = sign(x), and S(y) is the inverse of v(x) = xp(x) (y is increasing, by the ellipticity condition). The 
energy E (F q (R) = — V R E) is thus linear in p: 

E(K) = J d D rp(r)g q (\r~R\). (6) 

^From eq.(|5|), the effective Green's function, Q q , is integrated from (x ^ in light of the non-overlap assumption) 

VxG q (x) = s(qG)a S(z)^j, (7) 
|x| 

with z = |qC|/ao| x | D_1 ■ The Hessian of Q q is 

Qq,a = s (qG)a ^r-Bij, B i:j = % - j ^ f e^e,-, (8) 

|X| l ~\~ 

with = Xi/|x|. B has D — 1 eigenvalues equal 1, and one that equals —(D — 1)/(1 + ft) < 0. So, 

AS ? (x) = s(qG)(D i) ao S( jS )| x |-i T HL ( 9 ) 

(p,(z) = zp'/p). Also, 

^•(R)= / d°rp(r)g q . l3 (\r-R\). (10) 



Since S depends nontrivially on its argument, the dependence of F 9 , Q q , and E on q can be nontrivial. Among other 
things, the location of equilibria may depend on q. Take p to have a center of symmetry, which must then be an 
equilibrium point; when is it stable? If p has cubic symmetry the question hinges on the sign AE at the center (taken 
at the origin), because from the symmetry 

E ij (0)=D- 1 AE(0)S ij . (11) 



If /t and p have fixed signs, then from eqs.(|6|)(|9|) AE has the fixed sign s — sign(qfiGp) at all R for which there is no 
overlap. If s is positive, AE > 0, and there might be a true minimum (while AE < precludes a stable equilibrium). 
In particular, in the cubic-symmetric case the center is then per force a stable equilibrium. By the same token, the 
charge distribution p can be suspended (stably against translations) in the field of q, exemplifying statement (ii) 
below eq.(^). Thus, e.g., in the flow, Born-Infeld, and volume-minimization problems, where fiG < 0, p has to have a 
sign opposite to that of q to be able to suspend q. If the cubic symmetry of p is slightly disturbed there remains an 
equilibrium point (from continuity), which then may depend on q if the center of symmetry is lost. 

When p has lower-than-cubic symmetry, the center might still be a stable equilibrium: if Q qAj does not have a 
definite sign, Eij might be positive definite upon integration by eq.(|l0|). For example, in D > 2, if p is any planar 
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(two-dimensional) distribution with a four-fold symmetry, angular integration of B gives a matrix whose eigenvalues 
are all 1, but two that equal {fi, — D + 2)/2(l + fi). We get a stable minimum at the center if fi(z) > D — 2, for all z 
contributing to the radial integration. This is straightforwardly generalized to other symmetries. Clearly, because of 
the linearity in p stability of a composite configuration can be assessed by adding up Eij of the components. So, any 
combination of configurations that each gives a stable equilibrium also does so (e.g. a combination of any number of 
concentric, cubic-symmetric p(r)). 

What happens when q is decreased away from the large-g limit? In the general case, equilibrium might conceivably 
be lost, as q becomes comparable with the the "holding" charges, when one of the eigenvalues of Eij changes sign. 
For cubic-symmetric p, however, where eq.(p"l"|) is valid for all q, the answer depends on the sign of AE(0). In the 
test-particle limit, q — > 0, we have E — > qip. Since the Hessian of ip cannot have a definite sign at the center, we 
must have d q = q~ 1 AE(0) — * A<p(0) — > in the limit. It is likely that stability is lost in the test-charge limit with 
d q remaining positive all the way down to q = 0. Then, the center remains a stable equilibrium for all finite q. If 
d q changes sign and approach from below, then the same configuration has a stable equilibrium in a theory with 
G — > — G. (The possibility that d q = for a finite strech of q may be rejected from analyticity. In the linear case 
d q = for all q.) 

The above treatment (for large qj p) is applicable whenever we have the solution for the field ip B of B alone. For 
instance, all the above expressions still hold if B is any spherically symmetric charge (not overlapping with p) , with q 
its total charge. Exact solutions can also be found for any body with 1-D symmetry (plane, cylindrical, etc.). Another 
example: in |B] I derive the exact field for a pair of opposite point charges in the theory with p,(x) oc x D ~ 2 ; so its 
equilibria can also be studied in the limit of weak "holding" sources. 



IV. SOME IMPOSSIBILITIES 

Even if stable equilibria are possible in the nonlinear case, some elementary feats are still impossible (a) The 
impossibility to suspend a point, test charge in a single field stems from the fact that E = qip, and so AijEij = 
away from sources. As A is positive definite, Eij cannot be positive definite as stability requires. (I assume here and 
below that not all the Eij vanish at the equilibrium point. More generally, it can be shown that the first, relevant 
(even) derivative that does not vanish cannot be positive definite.) (b) It is impossible to balance stably a test point 
charge by two ip fields in a region where one of them has a constant gradient Vip' = a (e.g. levitate the charge 
in a ip field against earth-surface gravity). Here E = q(a ■ r + (p), and also satisfies Aij(Vip)Eij = away from 
sources, and so E cannot have a minimum, (c) Even if extended, test bodies can, in general be suspended, it is 
not possible to suspend a test dipole in a ip field. The force on a dipole p is (p • V)V</?. The translational energy 
is then E = p ■ Vip, and Aij(Vp)Eij = (p • V)(Aijp 4 . j ) — ip :i:i a ^ (p • V)ip, k . The first term vanishes everywhere 

away from sources; the second term vanishes at equilibria, because (p • V)V<£> is the force. Again, E^ cannot be 
positive definite at an equilibrium-no stable minimum, not even against translations, (d) Consider now a small test 
body that has both a charge q and a dipole p in equilibrium at a point where Vip ^ 0. Here E = qip + p • Vip. As 
in case (c) we have AijEij = —p ri j 9ip * 3 k (p ■ V)ip :k , but now the vanishing of the force implies (p • V)</? fc = —qip ih , 

so AijEij = qp,i,j p , k - Take, specifically, the orientation for which the moment on the dipole vanishes (also a 

requirement of equilibrium): p = —pSJip/\Vip\ (p > 0) calculated at the equilibrium point, and then consider stability 
against pure translations. For this value of p one finds that at the equilibrium point AijEij — p~ 1 q 2 \\7ip\ 2 ■ 
(putting ao = 1.) A stable equilibrium might then be possible only where fi is increasing. This is precluded, e.g., in 
the flow, and the volumc-minimizaton, problem, and in the modified dynamics, where /t is always decreasing, but not 
in Born-Infcld electrostatics. As we have only checked stability to translations, increasing fi does not guarantee full 
stability. 



V. SUSPENSION OF ALIGNED DIPOLES 

Neutral, dipolar bodies can be levitated stably, even in the linear case, provided the (constant-magnitude) dipole 
is forced to remain aligned with Vy. For example, under conditions of adiabaticity, precession of a spin (aligned 
with the dipole) about insures alignment. This fact is used, e.g., in the construction of magnetic traps for 
atoms and neutrons ||, and for macroscopic, magnetized tops Jic]]. Take then p = — pe, where e = Vtp/\Vp\ 
(restricted to equilibria at points where Vip ^ 0), and consider levitation of the body against a constant-force field 
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f. The total force is f + pV\Vip\, the energy (for translation plus realignment) is E(R) = — f • R p|V<^|, and 
Eij = — p\Vip\ ((p k,j — <P,*,i<P,*,j + l^Plv^.i,})) with z the direction of e (and is not summed over). As before, 
we require, as a necessary condition for stability, AijEij > at the equilibrium, where we can put (p^ A = — p /». 
After some algebra we get 

AijEij = -p-'lV^V^,^ - |f| 2 

+ p(f!-\t\ 2 )-p'\Vp\f z }. (12) 

In the linear case (p — 0) the expression in parentheses is positive (unless all the tp :i:j = 0.) and we restore p < 
as the necessary condition. In the case of suspension in a pure cp field (f = 0) it remains so in the nonlinear case. 
When f ^ 0, stability depends also on p and p! . Take, e.g., the case where f is along a symmetry axis of the field 
so fl = |f| 2 . Then, Ai 3 E i3 = —p~ 1 \Vip\~ 1 [p 2 ip t ^ k ip_ i:k - |f| 2 (l + p'\Vip\)]. With large enough p! it may be possible 
to get stable suspension even with p > 0. For example, when can such a dipole be levitated in the ip field of a point 
charge against a constant f? I find that this can be done if and only if one takes p > (dipole aligned with —Vip), 
and p satisfies at the equilibrium position /t^jjjf > 1 + (1 + p)/(D — 1). Of the problems listed above, the condition 
on p can be met only in the Born-Infeld theory; in the others pp' < 0. 
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